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^n : ABSTRACT 

in , 

r — [ ! We study the evolution of Super Star Cluster (SSC) winds driven by stellar winds 

^D I and supernova (SN) explosions. Time-dependent rates at which mass and energy are 

deposited into the cluster volume, as well as the time-dependent chemical composition of 
the re-inserted gas, are obtained from the population synthesis code Starburst99. These 
results are used as input for a semi-analytic code which determines the hydrodynamic 
properties of the cluster wind as a function of cluster age. Two types of winds are 
detected in the calculations. For the quasi-adiabatic solution, all of the inserted gas 
leaves the cluster in the form of a stationary wind. For the bimodal solution, some 
of the inserted gas becomes thermally unstable and forms dense warm clumps which 
accumulate inside the cluster. We calculate the evolution of the wind velocity and energy 
flux and integrate the amount of accumulated mass for clusters of different mass, radius 
and initial metallicity. We consider also conditions with low heating efficiency of the 
re-inserted gas or mass loading of the hot thermalized plasma with the gas left over from 
star formation. We find that the bimodal regime and the related mass accumulation 
occur if at least one of the two conditions above is fulfilled. 

Subject headings: Galaxies: star clusters — ISM: bubbles — ISM: HII regions — ISM 



Introduction 



Super star clusters are young compact objects observed in many starburst and interacting 



galaxies in a variety of wavelengths (IHgltzman et al. 



1992; 



Whitmore et al. 



2007 



1993; 



Galliano et al. 



Q'Connell et al. 



20081 ; 



Whelan et al 



1995 ; IMcIo et al.l[2005l ; ISmith et allboodlCilbert fc Grahaml 

201 ll ). With masses 10^ — 10^ Mq and ages < 10^ yr they are expected to include large numbers 



of massive stars which lose substantial fractions of their mass via stellar winds and supernova 
explosions. 



Chevalier Sz Cleggl (jl985l . hereafter CC85) studied the hydrodynamics of the gas re-inserted 
by massive stars into the cluster interior using an adiabatic spherically-symmetric model. They 
assumed that the mechanical energy of stellar winds and supernovae ejecta is thermalized in random 
collisions and the gas within the cluster is heated up to ~ 10^ K. The resulting high pressure drives 
the cluster wind for which CC85 found a stationary hydrodynamic solution. They assumed that 
the mass and the thermal energy are inserted uniformly at rates Msc and Lgc, respectively, into a 
sphere (cluster) of radius -Rsc- They showed, that under such assumptions, a stationary wind can 
only be obtained if the flow velocity equals zero at the cluster center and reaches the sound speed 
exactly at the cluster border. Super star clus ter winds we r e stu di ed further using a n alytical and 



numer i cal models by many autho rs including ICanto et al.l (|200d ) ; iRaga et al.l (|200ll ) ; ISilich et al, 
(|2003l ); iTenorio-Tade et al.l (|20od l. 



It was found by ISilich et al.l (120041 ) that the adiabatic approximation becomes inadequate for 
very massive and compact clusters. The authors showed that the stationary solution of the cluster 
wind does not exist for clusters with Lsc larger than a critical value Lcrit- This is because the total 
energy input rate, -^SC; is proportional to the cluster stellar mass, M^, while the energy losses from 
the hot gas due to radiation are proportional to M^ (since cooli r ig is p roportional to the second 



power of the gas density which is p roportional to M..). ISilich et al.l (|2004l ) showed how Lent depends 



on the star cluster parameters and lWiinsch et al.l (J2007l ) founded an approximate analytical formula 

for Lcrit- 



Clusters with Lsc > ^crit were studied by means of ID hydrodynamic simulations by lTenorio-Tagle et al. 



(|2007l ). who showed that such clusters evolve in the bimodal hydrodynamic regime. In such a case, 
the cluster is divided by the stagnation radius, Rst, into two qualitatively different regions. The 
stationary wind solution still exists in the outer region r > Rst, with the wind velocity being zero 
at Rst and reaching the sound speed at -Rsc- In the region r < Rst, on the other hand, the thermal 
instability sets in and random parcels of gas cool down to ^ 10^ K (further cooling is prevented by 
the intense stellar radiation). Consequently, the warm regions are compressed into dense clumps 
by repressurizing s hocks driven by the surrounding hot gas. Clusters in the bimodal regime were 
studied further by IWiinsch et al.l ()2008l ) who used 2D hydrodynamics to follow the clump forma- 
tion, and to estimate the fraction of the re-inserted matter which leaves the cluster as a wind and 
the fractio n which accumulates insid e the stagnation radius and possibly leads to secondary star 
formation (JTenorio-Tagle et al.ll2005l ). 



It was suggested that two-component supersonic recombina t ion line profiles of ten detected in 



young and massive SSCs (JGilbert &: Graham 



festation for such bimodal regime (JTenorio-Tagle et al 



2007 



dense HII regions overlapping vou ng SSCs (ISmith et al.ll200€ ) may present the obse rvational mani 



201C 



Becd I2OO8I : iHenrv et al.1 120071 ) and compact 



Silichet al.ll2007l . l2009l ). In both cases 



the calculations require the shocked gas temperature to be lower than that predicted by the CC85 
model as it is al so the case when the mo del predicted diffuse X-ray emission is compared to the 
observed values (jStevens &: Hartwelll l2003l ) . Two different processes which may decrease the inter- 
cluster gas temperature have been discussed in the literature: the efficiency with which the kinetic 
energy of stellar w inds and SNe is thermalized, and the additional mass loading i n to the hot gas 



insid e the cluster dSt evens fc Hartwelll l2003l : iMelioli fc de Gouveia Dal Find |200J; IWiinsch et al. 



2007 



Silich et alj 1200 7. 2009, l20ld ). In this work we do not discuss details related to those two 



processes, however, we introduce two free parameters r^he and r^mi and show how the results depend 
on their values. 

Previous works on clusters in the bimodal regime use the energy and mass deposition rates 
Lsc and Msc as free parameters. In this work, we calc ulate time-dependent Lsc(^) and Msc(t) 
using the stellar population synthesis code Starburst99 ([Leitherer et al.lll999l ) for a cluster with 
a given stellar mass, M*, and initial stellar metallicity, Zq. Subsequently, we insert Lscit) and 
Mscit) into our semi-analytic code to determine the evolutionary properties of the cluster wind. 
We also calculate whether the cluster spends some time in the bimodal regime and estimate the 
amount of re-inserted gas which becomes thermally unstable and accumulates inside the cluster. 
The Starburst99 code also provides us with the time evolution of the re-inserted gas chemical 
composition. The chemical composition is an important parameter as the cooling rate depends on 
it. This work effectively replaces the three functions of time Lgc(i), Mscit) and Z{t) (metallicity 
of the cluster wind), with the two constant parameters: mass of the star cluster M^ and its initial 
metallicity Zq. 

The paper is organized as follows: in S|2] we describe the semi-analytic code used for the 
calculation of the cluster wind and the way how it utilizes results of the Starburst99 code. In 
^ we show results for a reference model with M^ = 10^ Mq and Rsc = 3 pc ( §3.ip and give 
the dependence of results on the cluster mass, the cluster radius and the initial stellar metallicity 
( §3.2p . In f|l]we summarize our conclusions. 



2. The cluster v^ind 



In this section we specify the assumptions used in the semi-analytic model of the cluster wind 
and formulate its basic equations. We also describe how the model equations are integrated and 
the properties of the bimodal solution (e.g. Rst, -^crit and M^ind) determined. Finally, we describe 
how the wind model utilizes the output from the Starburst99 code. 
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Fig. 1. — Energy (solid line, left y-axis) and mass (dashed line, right y-axis) deposition rates 
calculated by the Starburst99 code for the reference model M^, = 10^ M©, i?sc = 3 pc, Zq =Lq., 
r]i,c = 1 and r/mi = 0. 




Fig. 2. — The evolution of the metalhcity of the hot shocked gas inside the cluster with Zq = Zq. 
The solid and dashed lines show the metallicity without (rj^i = 0) and with (rj^ni = 19) mass 
loading, respectively. 



2.1. Assumptions and basic equations 



We consider a spherical cluster of radius Rsc whose stars deposit mass and energy at rates 
Msc and Lsc, respectively. We assume, similar to CC85, that mutual collisions of stellar winds 
and SNe ejecta, and collisions with gas left over from the formation of the first stellar generation 
and with gas returned by pre- main sequence stars via outflows, result in the production of hot 
gas which occupies most of the cluster volume. Therefore, we model these processes by inserting 
mass and energy uniformly into the whole cluster volume with deposition rate densities qm and qe. 
In order to account for the uncertainties related to the thermalization of the mechanical energy 
of the inserted gas we introduce the heating efficiency, r/he & (Ojl); denoting the fraction of the 
mechanical energy of stellar winds and SN ejecta which is converted into the thermal energy of the 
hot gas. Furthermore, we assume that a substantial amount of gas was left over after the formation 
of the first generation of stars and that this gas may evaporate and be dispersed into the hot gas. 
Another contribution to the mass of hot gas inside the cluster comes from outflows of pre-main 
sequence stars which are not included in the Starburst99 code. Indeed, the mass lef t over from star 



formation and the T Tauri multiple outflows, such as the jet from RW Aurigae (JBacciotti et al. 



19961 ). should make a substantial contribution to the mass available for mass loading. We describe 
these processes by the mass loading factor, rymi G (0,oo), which gives the amount of the loaded 
mass relative to Msc. The total mass injection rate then is Min = Msc + ^miMsc = (1 + ^mi)Msc- 
It is assumed that the metallicity of the loaded gas is the same as the initial stellar metallicity, Zq. 



The spherically symmetric hydrodynamic equations describing the wind flow are (CC85. ISilich et al. 
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where 7 is the adiabatic index and p, u and P are wind density, velocity and pressure, respectively. 
Mass and energy deposition rate densities qm and qe are 



Qm 
Qe 



3(l+r;^l)Msc 



4nRfc 






(4) 



for r < Rsc and qm = qe = for r > Rsc- The energy equation (^ includes the cooling term 
Q = nineA{T, Z) where rii = n^ = p/ Pi are the ion and electron numb er densities, Z is the gas 
metallicity and A(T, Z) is a cooling function tabulated bv iPlewal (|l995l ). We use pi = 14/llm,H 
neglecting the contribution of heavy elements. 



Several interesting properties may be derived directly from equations ([I]) - ([3|) (see lSilich et al 



2004 . for details). Firstly, the stationary solution exists only if the wind velocity, u, reaches the 
sound speed exactly at the cluster border. Secondly, a relation between the temperature Tgt and 
the density pst at the stagnation radius can be derived 



Pst 



qe - qmclt/il 



1) 



KTsuZ) 



1/2 



(5) 



where Cgt is the sound speed at the stagnation radius. Furthermore, it has been shown bv lTenorio-Tagle et al 
(|2007l ) that if the cluster is in the bimodal regime, i.e. if R^t > 0, the pressure at the stagnation 
radius Pst = {kpstTst)/tJ-a reaches the maximum value Pmax = max(Pst(Tst)), where k denotes the 
Boltzmann constant and ^a = 14/23m,H is the mean mass per particle. 



2.2. Integration procedure 

The wind solution is found by the following procedure. At first, it is assumed that R^t = and 
an attempt to find Tgt is made. Equations ([I]) - ([3]) are repeatedly numerically integrated from r = 
to Rsc with Tst varying in the interval (0, T^) where T^ = (7 — l)fiaQe/i'ykqm) is the adiabatic wind 
central temperature. The central density pst is calculated from equation ([5]). Then, the bisection 
method is used to find Tst for which the sonic radius Rson (defined as n(i?son) = Cs(-Rson)) is equal 
to itse- 
lf this attempt fails (i.e. no initial conditions at r = for which -Rson = Rsc exist), it implies 
that Rgt > and the cluster is in the bimodal regime. In such a case , the value of Tst is define d 
by the condition that the function Pst(?st) has its maximum Pmax (ITenorio-Tagle et al.l 120071 ). 
Therefore, the temperature at the stagnation radius is found using the golden section method and 
it is used as the initial condition for integrating equations ([I]) - ([3]). Then, similarly as in the 
previous case, Rst is varied and the bisection method is used to find the solution which satisfies the 
condition Rson = Rsc- 

Once all the initial conditions {Rst, Pst and Tst) are known, radial profiles of the wind density 
p, velocity u and temperature T can be obtained by integrating equations ([I]) - ([3]) in the interval 
(i?st) 10 Rsc- The semi-analytic model is unable to de scribe the inner ther mally unstable region 
with r < Rst- However, 2D hydrodynamic simulations (jWiinsch et al.ll2008l ) have shown that the 
temperature and the density of the hot gas in this region are close to uniform and stay constant 
with time. The deposition of mass into this region is balanced by the formation of dense warm 
clumps which tend to accumulate in this region. Therefore, we assume that the hot gas in the 
central region r < i?st has zero velocity, uniform density pst and temperature Tjt, and that all gas 
inserted into this region accumulates there. Finally, the critical luminosity Lent is determined by 
repeating the above procedure and searching for the lowest mechanical luminosity Lgc for which 
Rst > 0. 



2.3. Starburst99 outputs used in the wind model 



The stellar population synthesis code Starburst99 (jLeitherer et al.lll999l ) calculates a set of 
stellar evolution models for a given population of stars and determines their collective properties. 
In this work, the total mass loss rate from stellar winds and SNe type II ejecta is used as the mass 
deposition rate, Msc and the total stellar wind and SNe ejecta power as the energy deposition 
rate, -Lsc- All Starburst99 simulations used in this work are set up with the following parame- 
ters: star form ation is instan taneous with the fixed stellar mass M^,; the standard Kroupa Initial 
Mass Function ( Kroupal 120011 ) with two power-laws {dN/dm ~ m~^'^ between 0.1 and 0.5 M0and 
dN/dm ~ m~^'^ between 0.5 and 100 Mq) is used; the supernova cut-off mass is equal to 8 Mq; 
stellar evolutionary t racks are Geneva with high mass loss; and the wind model is evolutionary (see 



Leitherer et al.lll992l for details). The evolution of Msc and Lgc for the reference model (see ^3.ip 
is shown in Figure [TJ We have followed the first 40 Myr of the cluster evolution. This period is 
long enough to cover the life time of all massive stars even in cases with initial stellar metallicities, 
Zq, different than Zq, discussed in ^3.2[ We do not consider here the period after the last massive 
star explodes (this moment is visible as a sudden drop of Msc and Lsc at 37 Myr in Figure [T]) . 

Starburst99 also provides the chemical composition of the re-inserted matter by specifying 
mass loss rates for: H, He, C, N, O, Mg, Si, S and Fe. Thus one can calculate the injection rate 
for seven elements heavier than H and He 

Fe 



Mi^etals = Y^ Mj (6) 



where Mj is the mass deposition rate of the j-th element. It is assumed that the injected gas is 
rapidly mixed with the mass loaded gas. The metallicity of the cluster wind, Z, used in equation 
©is 

(1 + VmlWsc 

The evolution of Z in the cluster with Zq = Zq for different values of T^mi is shown in Figure [2j 
Taken together, the model utilizing Starburst99 results includes five parameters: M^, Rsc, Zq, -q^e 
and r/mi. The semi-analytic wind model on its own includes six parameters: Mgc, Lsc, Z, Rsc, 
r/he and r/mi- Here we assume that the first three of them (Mgc, Lsc, Z) are functions of the star 
cluster age. We keep the heating efficiency, r/he, and the mass loading coefficient, rjrai, constant, 
despite they may change with time as the number of massive stars and the amount of gas left over 
from star formation decrease. 



3. Results 

In the first part of this section ( ^3.ip . we describe in detail results for our reference model 
whose parameters are chosen to represent a typical SSC. Since the heating efficiency, ??he) and the 



mass loading factor, rj^i, are free parameters, we show results for three different combinations of 
them. In section ^3.21 we show how the most important results (the existence of the bimodal regime 
and the amount of the accumulated mass) depend on the cluster mass, the cluster radius and the 
initial metallicity of stars and the mass- loaded gas Zq. 



3.1. The reference model 

We calculate the evolution of a wind driven by a cluster with a stellar mass M^ = 10^ Mq, 
radius i?sc = 3 pc and initial stellar metallicity Zq = Zq = 0.02. We explore three combinations of 
r/he and ??mi (see Tabled]). In model A, ther e is no mass l oading and the heating efficiency is 100%. 
Model B is chosen to be in agreement with lSilich et al.l (J2007l . |2009J) who have obtained r]^c — 5%, 
in order to fit the parameters of the compact HII regions observed around 11 SSCs selected in the 
central zone of M82. In model C, the mass loading factor, r]^\ = 19, is set to give the same value 
of ^ 00 5 as in model B, where 



Vr 



r],oo 



2?/ho-^sc 



1/2 



(8) 



(1 + ??ml)Msc 

is the adiabatic wind terminal speed corrected for effects of heating efficiency and mass loading. 

Figure [3] compares the time evolution of the critical luminosity, Lent, with the star cluster 
mechanical luminosity, Lsc- In model A the star cluster mechanical luminosity is always below 
the critical value, Lsc < -^criti and thus all gas re-inserted by stars leaves the cluster as a wind. 
On the other hand, models B and C present periods with Lgc > -^crit when clusters evolve in the 
bimodal regime. The beginning and the end of these periods are shown in Table [1] in columns tbs 
and tbe! respectively. Even though models B and C have the same Fr?,ooj the period of bimodality 
is longer in model C. This is because, due to mass loading, the density of the thermalized plasma is 
larger in model C and it results in a higher cooling rate that favors thermal instabilities and mass 
accumulation. 

The stagnation radius, Rst, for the three models is shown in FigureHl For model A, it is always 
at the cluster center, while in the other two cases, Rst reaches a substantial fraction of i?sc when 
the clusters evolve in the bimodal regime. This implies that the amount of mass accumulated in the 



Model 


Vhe 


^ml 


ibs [Myr] 


the [Myr] 


Mace [Mq] 


Min [Mq] 


A 
B 
C 


1 

0.05 
1 




19 


2.4 
1.3 


11.1 
16.9 



5.8 X 10^ 
1.8 X 10^ 


1.8 X 10^ 
1.8 X 10^ 
3.7 X 10^ 



Table 1: The reference model calculated with different ryhe and 77^1. Columns 4 and 5 denote the 
beginning and the end of the period of bimodality (L^c > ^crit)- Columns 6 and 7 show the amount 
of mass accumulated inside the cluster. Mace, and the total amount of mass, Mjn, supplied into the 
cluster by stars and mass loading, respectively. 
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Fig. 3. — The evolution of the critical luminosity, Lcrit, for models A (dashed), B (dotted) and 
C (dash-dotted). The Lcrit curves are compared to the cluster mechanical luminosity Lsc (solid). 
Periods during which the cluster evolves in the bimodal regime are denoted by arrows. 
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central zones of the cluster may be significant if the heating efficiency is low or the mass loading is 
large. It is because the mass accumulation rate is Ma.cc = -^in(-Rst/-Rsc)^ where Min = (l+??mi)-^sc 
is the rate at which mass is supplied into the cluster by stars and mass loading. For example, the 
amount of the accumulated matter, M^cc = /j ^° -^accdi, is about one third of the total mass 
supplied into the cluster, Min = /j'"' -^indt, in the case of model B and about one half of Min in 
the case of model C (see Table [1|) . 

Note that strong radiative cooling also affects the star cluster wind mechanical output rate, 
-^wind = 47r/9ur^(u^/2 + H), where H is the enthalpy. Figure [5] shows that in the bimodal regime 
it falls well below the star cluster mechanical luminosity, Lsc (model C), and below the heating 
efficiency reduced star cluster mechanical luminosity, r/he-^sc (model B). This implies that the 
"true" energy output and thus the impact of SSCs on the ambient ISM may be much smaller 
than one would expect from star cluster synthetic models like Starburst99. Note also that the 
star cluster wind terminal speed is in such cases smaller than that expected from the star cluster 
synthetic models, see Figure [H It compares the wind terminal speed Voo (measured from semi- 
analytic models at r = 10 -Rsc) to the heating efficiency and mass loading corrected adiabatic wind 
terminal speed V^^oo- The difference between the two, significant mainly during bimodality periods, 
is due to the radiative energy losses from the wind. 



3.2. Dependence on the stellar cluster parameters 

In this section we explore how our results depend on the cluster parameters running models A, 
B and C for clusters with different masses {M^ = 10^,10^ and 10^ Mq), different radii (-Rsc = Ij 
3, 10 and 30 pc) and taking into consideration the variation of the re-inserted gas metallicity. The 
results of the calculations for stellar clusters with different masses and radii, when the re-inserted 
and the ablated gas metallicity were fixed to the solar value are presented in Figure [71 This 
figure compares the calculated critical mechanical luminosities, Lcrit, to the star cluster mechanical 
luminosity obtained from the Starburst99 synthetic model. Figure \7\ shows that clusters with 
r/he = 1 and r/mi = never evolve in the bimodal regime. On the other hand, models with low 
heating efficiency or large mass loading exhibit periods of bimodality (see Tables [2] and [3]) . In 
the extreme cases the amount of mass accumulated inside the cluster. Mace, may reach 70% of 
the re-inserted and ablated mass, as it is the case when rj^i = 19 and M^ = 10^ Mq. Note that 



Wiinsch et al.l (|2007l ) derived an approximate analytic formula for Lent which predicts that Lcrit 
is in direct proportion to the size of the cluster, .Rsc- This is in excellent agreement with our 
semi-analytic results. Note also that both Lsc and Mgc are linearly proportional to Mj, resulting 
in Lcrit independent of M+. Thus, Lent defines the critical cluster mass. Merit, and clusters evolve in 
the bimodal regime if AL > Merit- This linear dependence may be broken if the cluster IMF varies 
with the cluster mass, or if more massive clusters are formed in a different more abrupt process 
compared to low mass clusters. However, in this paper we explore consequences of an abrupt cluster 
formation with a given IMF. Discussion of their dependence on the cluster mass exceeds the scope 
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Fig. 4. — The evolution of the stagnation radius, Rst, for models B (dotted) and C (dash-dotted). 
The stagnation radius is always zero in model A. 
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Fig. 5. — The evolution of the wind mechanical output rate, -Lwind) for models A, B and C is shown 
by thick dashed, dotted and dash-dotted lines, respectively. The thin solid curve shows the star 
cluster mechanical luminosity, L^c- It overlaps with the thick dashed curve, because for model A, 
cooling from the hot wind is negligible and L^ind = -^sc- The thin double-dashed line represents 
the heating efficiency reduced energy deposition rate, r/hc-^sC) for niodel B with rji^^ = 0.05. 
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Fig. 6. — The evolution of wind terminal speed, Foo- Thick dashed, dotted and dash-dotted lines 
show the calculated wind terminal speed in the case of models A, B and C, respectively. Thin lines 
display Fr;,oo (see equation [8]) , for model A (solid) and models B and C (double dashed). 
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of this paper. The results of the calculations for clusters with different masses and radii in the 
case when the inserted gas metallicity is solar are summarized in Tables [2] and [3l The tables show 
that even in the case of low heating efficiency or lar ge mass load i ng, cl usters evolve in the bimodal 



regime only for some time, as it was suggested in ISilich et alj (|2009l ). The length of the period 
of bimodality and the amount of accumulated mass are larger for clusters with smaller radii and 
larger masses. 

Another parameter which may affect properties of the star cluster driven outflows is the re- 
inserted gas metallicity. In the case of instantaneous star formation, the metallicity of the re-inserted 
matter changes a lot, as it is shown in Figure [21 T his should change the coo ling rate and thus the 



critical mechanical luminosity, Lcrit, significantly (jTenorio-Tagle et al.ll2005l ). In order to explore 
how our results depend on this parameter, we have varied the initial stellar and the loaded gas 
metallicity, Zq, in our reference models A, B and C. Three different values of Zq were used for the 
calculations: Zq = 0.05Zq, Zq = Zq and Zq = 2.0Zq. The top left panel in Figure [8] shows the 
trends of the wind metallicity, Z, calculated from equation[71 In all cases without mass loading (solid 
lines in Figure [8]) the metallicity of the thermalized plasma grows rapidly after the first supernova 
explodes, reaches about 10 times the solar value, and then decreases gradually reaching about 3 
times the solar value after ~ 20 Myr. In the case with mass loading, the maximum metallicity 
never reaches 10 times the solar value. This is because in this case the re-inserted matter mixes 
continuously with a large amount of the ablated gas. The calculated critical luminosities, -Lcrit, 
are then compared with the star cluster mechanical luminosities, Lsc (top right, bottom left and 
bottom right panels in Figure [8] for cases A, B and C, respectively). Models without mass loading 
and r/he = 1 never enter the bimodal regime (see top right panel). Note that relative abundances 
of species in the re-inserted matter differs from solar values. This implies that the cooling function 
using scaled solar composition Z may give somewhat different cooling rates that that calculated 
from individual species separately. This, however, does not change our results significantly, since 
the main coolants (C and O) are also dominant ingredients of Z. 



4. Conclusions 

We used our semi-analytic spherically-symmetric code together with the stellar population 
synthesis model Starburst99 to study the time evolution of Super Star Cluster winds. 

Two physical processes which could affect the hydrodynamics of the star cluster winds signif- 
icantly and cannot be studied in the semi-analytic approach in details, the heating efficiency and 
mass loading, are parameterized with two constant parameters r/he and r/mi- We also search how 
our major results depend on the metallicity of the re-inserted matter. 

The calculations show that strong radiative cooling becomes a crucial issue when the wind is 
mass loaded or the thermalization efficiency (and thus the fraction of the star cluster mechanical 
luminosity which drives the outflow) is small. In these cases (our reference models C and B, 
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1.8 X 10^ 


10^ 


1 





- 


- 





1.8 X 10^ 


10^ 


1 





- 


- 





1.8 X 10^ 


10^ 


0.05 





3.1 


5.1 


1.1 X 10^ 


1.8 X 10^ 


10^ 


0.05 





2.4 


11.1 


5.8 X 10^ 


1.8 X 10^ 


10^ 


0.05 





1.6 


17.3 


1.0 X 10^ 


1.8 X 10^ 


10^ 


1 


19 


2.1 


9.8 


8.0 X 10^ 


3.7 X 10^ 


10^ 


1 


19 


1.2 


16.9 


1.8 X 10^ 


3.7 X 10^ 


10^ 


1 


19 


0.0 


36.8 


2.8 X 10^ 


3.7 X 10^ 



Table 2: Clusters with different stellar mass, M^, heating efficiency, r/hej and mass loading rj^i. The 
cluster radius is i?sc = 3 pc for all these models. Columns 4-7 have the same meaning as in 
Table H 



^sc [pc] 


Vhe 


^ml 


ibs [Myr] 


the [Myr] 


Mace [Mq] 


Mi„ [Mq] 


1 


1 





- 


- 





1.8 X 10^ 


3 


1 





- 


- 





1.8 X 10^ 


10 


1 





- 


- 





1.8 X 10^ 


30 


1 





- 


- 





1.8 X 10^ 


1 


0.05 





2.1 


13.4 


8.0 X 10"^ 


1.8 X 10^ 


3 


0.05 





2.4 


11.1 


5.8 X 10^ 


1.8 X 10^ 


10 


0.05 





2.8 


9.2 


3.1 X 10^ 


1.8 X 10^ 


30 


0.05 





3.1 


5.1 


1.1 X 10^ 


1.8 X 10^ 


1 


1 


19 


0.2 


24.2 


2.3 X 10*^ 


3.7 X 10^ 


3 


1 


19 


1.2 


16.9 


1.8 X 10*^ 


3.7 X 10^ 


10 


1 


19 


1.8 


12.4 


1.3 X 10^ 


3.7 X 10^ 


30 


1 


19 


2.1 


9.8 


8.0 X 10^ 


3.7 X 106 



Table 3: Clusters with different radius, Rsc, heating efficiency, 7?he) and mass loading r^mi- The 
cluster stellar mass is M^ = 10^ Mq for all these models. Columns 4-7 have the same meaning 
as in Table [TJ 
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Fig. 7. — Thin lines show the evolution of the critical luminosity, Lcrit, for models A, B and C (the 
lines have the same meaning as in Figure [3]) and different cluster radii, Rsc^ marked in the figure. 
Thick solid lines display the mechanical luminosity, -Z^sC) foi' clusters with different stellar masses 
as denoted in the figure. 
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Fig. 8. — Dependence on the stellar metallicity and the metallicity of the mass loaded gas Zq. Top 
left panel shows the metallicity of the supplied gas, Z, for models with rj^i = (solid, models A 
and B) and ?7mi = 19 (dashed, model C). Other panels compare the evolution of Lsc and Lent with 
different Zq for models A (top right), models B (bottom left) and models C (bottom right). In all 
panels, colors of curves represent Zq: Zq = 0.05 Zq (red), Zq = Zq (green) and Zq = 2.0 Zq (blue). 
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respectively) the evolutionary tracks of the star cluster winds show periods of bimodality. During 
these periods only some fraction of the re-inserted and loaded gas leaves the cluster as a wind. The 
rest of the re-inserted matter cools down rapidly, becomes thermally unstable and is accumulated 
in the central region of the cluster. The duration of these periods depends on the star cluster 
parameters rj^c and 77^1. Periods of bimodality are longer in the case of more massive clusters with 
smaller radii. However, they become progressively shorter as the mass loading drops or the heating 
efficiency grows. The bimodal regime vanishes in the cases when heating efficiency is large and mass 
loading is insignificant. In the simulations which include mass loading, the stellar metallicity does 
not affect significantly neither the duration of the bimodal regime nor the amount of re-inserted 
mass which accumulates inside the cluster. Models with low heating efficiency are more sensitive 
to the metallicity of the re-inserted matter. 

We conclude that the second stellar generation may be formed in massive and compact stellar 
clusters from thermally unstable parts of stellar winds and the mass loaded gas in their central 
parts. Low heating efficiency ?7hc = 0.05 leads to the second stellar generation heavily enriched 
with He-burning products. However, its total mass is a few percent of the first generation only. 
High value of mass loading r/mi = 19 results in the massive second stellar generation, however, its 
metallicity is only slightly higher than that of the first generation. 
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Zo [Zo] 


??he 


^ml 


tbs [Myr] 


the [Myr] 


Mace [Mq] 


Min [Mq] 


0.05 


1 





- 


- 





1.8 X 10^ 


1.0 


1 





- 


- 





1.8 X 10^ 


2.0 


1 





- 


- 





1.8 X 10^ 


0.05 


0.05 





3.1 


12.9 


7.5 X 10^ 


1.8 X 10^ 


1.0 


0.05 





2.4 


11.1 


5.8 X 10^ 


1.8 X 10^ 


2.0 


0.05 





1.9 


9.3 


3.7 X 10^ 


1.8 X 10^ 


0.05 


1 


19 


2.8 


15.0 


1.6 X 10^ 


3.6 X 10^ 


1.0 


1 


19 


1.2 


16.9 


1.8 X 10^ 


3.7 X 10^ 


2.0 


1 


19 


0.4 


17.0 


1.8 X 10^ 


3.6 X 10^ 



Table 4: Models with different stellar metallicity, Zq, heating efficiency, r/he, and mass loading j^mi- 
Other cluster parameters are i?sc = 3 pc and M^ = 10^ M©. Columns 4-7 have the same meaning 
as in Table [TJ 



